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ABSTRACT 


The  usual  mathematical  formulation  of  availability- 
assumes  an  exponential  distribution  for  failure  and  repair 
times.   While  this  is  sometimes  correct  for  reliability,  it 
is  not  likely  to  be  for  maintainability.   This  study  was 
conducted  to  verify  that  the  lognormal  and  gamma  distribu- 
tions are  suitable  descriptors  for  corrective  maintenance 
repair  times,  and  to  estimate  the  difference  caused  in 
assuming  an  exponential  distribution  for  availability  and 
maintainability  calculations  when  in  fact  the  distribution 
is  lognormal.   Forty- six  sets  of  data  of  electronic  and 
mechanical  systems  and  equipments  were  analyzed  using  the 
methods  of  probability  plotting  and  statistical  testing  for 
distributional  assumptions. 


TABLE    OF    CONTENTS 


I. 


II 


INTRODUCTION 


OVERVIEW 


B.   PURPOSE  AND  APPROACH 


1.  Objectives  -------------  9 

2.  Systems  and  data  analyzed   -----  13 

3.  Analysis  approach   ---------  10 

STATISTICAL  CONSIDERATIONS  ---------  14 


A.   ROUGH  CHECK  METHODS 


1.  Boxplot   -  -  -  -  -  -  . 

2.  Histogram   ----- 

B.   STATISTICAL  DISTRIBUTIONS 


1.  Exponential  Distribution 

2.  Lognormal  Distribution 

3.  Gamma  Distribution   -  - 


C.   PROBABILITY  PLOTTING    -  -  - 

1.  Lognormal  Distribution 

2.  Exponential  Distribution 


3.   Gamma  Distribution 


D.   STATISTICAL  TESTS  FOR  DISTRIBUTIONAL 


ASSUMPTIONS 


14 
14 
16 
17 
17 
18 
19 
20 
21 
22 
22 

23 


1.      The   Chi-Sauared   Goodr.r ss-Of -Fi- 


Test 


2U 


2.  The    Kclmogorov-Smirno*    Test      -   -    -    -  26 

3.  The    W    Test    for   Lcgnormal 

Assumption   -------------  28 

III.       DATA    ANALYSIS    RESULTS       -----------  30 

A.  RESULTS    OF    TESTS    FOR    DISTRIBUTIONAL 
ASSUMPTIONS       ------------  30 

1.  Summary  of   results    ---------  30 

2.  Discussion   of   specific:    cases    -   -    -    -  39 

B.  ERRORS    IN    CALCULATED    PARAMETERS       -    -    -    -  114 
1.      Error    in    HTTP    when   assuming    an 

exponential    distribution 

instead  of   lcgnormal    --------  44 

C.  CONCLUSIONS       --------------  46 

APPENDIX  A.   -------------------  47 

APPENDIX  3.   -------------------  94 

APPENDIX  C   Computer  Programs  for 

Data  Analysis  ------------  99 

LIST  OF  REFERENCES  ----------------  106 

INITIAL  DISTRIBUTION  LIST   ----------  109 


LIST  OF  TABLES 

1.  System  /  Equipments  analyzed.  --------  11 

2.  Total  data  analysis  results  ---------  31 

3.  Lognormal  and  Exponential  results  ------  32 

1*.   Gamma  distribution  tes+  results  -------  33 

5.  Percentage  error  in  HTTR  when  assuming 
an  exponential  instsad  cf  lognDrmal 
distribution   ----------------  45 

6.  To-'-.al  results  ----------------  94 

7.  Pump  data  results  --------------  95 

8.  Report  data  results  -------------  96 

9.  Field  data  results   -------------  97 

10.   Total  percent  results  ------------  98 


ACKNOWLEDGMENT 


The  author  wishes  to  express  his  appreciation  to  the 
individuals  who  have  given  their  time  and  support  in  this 
effort.   Specifically,  I  wDuld  liks  to    recognize  Professor 
H.  3.  Kline  for  his  guidance,  and  assistance  and  patience; 
Professor  D.  P.  Gaver  as  the  second  reader;  the  members  of 
the  Electricite  de  France  ,  Rome  Air  Development  Center  and 
United  Kingdom  Ministry  of  Defence  who  provided  us  data;  and 
my  wife  Sebnem,  for  her  help  and  patience  during  this  study. 


I. INTRODUCTION 

A.  OVERVIEW 

The  effectiveness  of  a  system  depends  not  only  on  its 
ability  to  meet  its  specified  performance  requirements,  but 
also  on  its  ability  to  perform  when  needed,  for  the  duration 
of  its  assigned  missions,  and  for  its  operational  lifetime. 
The  technical  disciplines  concerned  with  these  time-related 
system  characteristics  are  reliability,   maintainability  and 
logistic  support. 

Commonly  used  methods  for  prediction  and  demonstration 
of  corrective  maintenance  repair  times  assume  the  validity 
of  the  lognormal  distribution. 

In  order  to  help  in  focusing  attention  on  these  matters, 
a  statistical  analysis  on  data  sets  of  demonstrated  and 
field  repair  times  has  been  conducted  as  part  of  a  study  on 
the  application  of  the  lognormal  and  gamma  distributions  to 
corrective  maintenance  data.   Ihe  results  of  the  analysis 
are  given  in  this  thesis.   A  previous  study  of  electronic 
systems  and  equipments  concluded  that  the  lognormal  distri- 
bution was  a  suitable  descriptor  for  repair  times  [Ref.  1]. 

B.  PURPOSE  AND  APPROACH 
1 ■   Ob  iectives 

The  objectives  of  this  study  have  been; 


(a)  To  test  the  exponential,  lognormai,  and  gamma 
distributions  as  descriptors  of  repair  time  data  on  mechani- 
cal and  electronic  systems  and  equipments, 

(b)  to  verify  that  the  lognormai  or  gamma  distri  - 
bution  is  a  suitable  descrip-or  for  repair  times, 

(c)  to  estimate  the  percentage  error  caused  in 
assuming  exponentiali ty  for  availability  and  maintainability 
calculations  when  in  fact  the  distribution  is  lognormal. 

2  •   Syste ms  and  Data  Analyzed 

Forthy-six  sets  of  repair  time  data  for  electronic 
and  mechanical  systems  /  equipments  vsre  analyzed  (Taole  1). 
Data  for  the  electronic  systems  (Sets  8-20)  cama  from  formal 
maintanability  demonstration  test  reports  furnished  to  us  by 
the  Rome  Air  Development  Center,  U.S.  Air  Force.   Mechanical 
equipments  included  field  repair  time  pumps  used  in  French 
nuclear  electric  power  stations  (Sets  1-7),  and  additional 
26  sets  of  field  repair  data  from  British  Aircraft,  provided 
by  the  United  Kingdom  Ministry  of  Defence,  representing 
elapsed  times  (Sets  21-33),  and  man  hours  (Sets  34-46).   >Io 
detailed  reports  are  available  for  the  mechanical  items, 
3-   Analysis  Approach 

Three  approaches  were  used  for  testing  data. 
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Table    1 


SYSTEMS       /       EQUIPMENTS    ANALTZED 


EM  N 


Set 

No 

— r 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

15a 

16 

17 

18 

19 

20 

21, 

,34 

22 

,35 

23 

,36 

24, 

,37 

25, 

,38 

26 

,39 

27 

,40 

28 

,41 

29, 

,42 

30, 

,43 

31  , 

,44 

32, 

,45 

33, 

,46 

M 


FfencH  F 

French   F 

Condense 

Heater   D 

Main   Fee 

Reactor 

Condense 

AN/GRT-2 

AN/ARC- 1 

AN/ARC-1 

AN/ALQ-1 

CGT    MAX. 

AN/TPN-2 

AN/MST-T 

AN/UYK-1 

Revised 

Integrat 

CCG-NBDT 

AN/GXS-2 

AN/ARN-1 

AN/FCC- 9 

Fuselage 

Flying   Z 

Landing 

Hydra uli 

Ice  and 

Power    PI 

Engine   S 

Mam    Rot 

Main    Rot 

Electric 

Co  mmur.ic 

Surveill 

Automat  i 


esse 
esse 
r  Ex 
rain 
dwat 
Cool 
r  Va 
0 

62  (7 
6  3  (V 
25    I 

MPU 
8  Ra 
1  El 
4C 

AN/U 
ed    M 

Com 
V 

01 
8 

Str 
oner 
Gear 
cs 

Rain 
ant 
tart 
or  H 
or  D 
al  P 
atio 
ance 
c  Fl 


nKeim    Pumps   D~own~Time 
nheiia    PumDs    Repair    Time 
traction    Pump 

Pump 
er  Punp 
ant  Pump 
cuum  Pump 

Radio  Set  Intermediate 
UHF  Modern  Radio  S?4: 
EREC  ,  Direction  Finder  Set 

and  RSS  Terminal  Orqanize 
dar  Bsacon  Transponder 
ectronic  Warfare' Training  Si 


YK-143 
icrowive  G/G 

pass  Ears 


Subs  ystem 


ffi 


Itiplrxer 

ucture 

ols 


Prot action 

ing 

ead  and  Blades 

rive 

ewer  Supply  and  Distribution 

ns 

and  Search 
ight  Control 


Note:  In  sets  21-45,  the  two  numbers  represent  two 
part  data-  The  first  part  is  ELAPSED  TIME  and  second  part 
is  MAN  HOURS. 
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(a)  Rough  check  using  a  technique  called  Boxplot. 
This  gives  a  rough  idea  that  the  data  comes 
from  symmetric  (normal)  or  skewed  distributions 
lognormal,  exponential  or  gamma). 

(b)  Plotting  the  data  using  appropriate  probability 
paper. 

(c)  Statistical  tests  of  significance  for  assumed 
distributions. 

The  following  procedure  was  used  for  analyzing  the 


dat  a, 


i.      Sketch   boxplot    and   histogram  of   data    and 
determine    if    data    appears   to   come   from    a    symmetric   or    skewed 
distribution.       Obtain    an    idea   about    center    and   dispersion    of 
dat  a. 

ii.      Plot    the    data   on    appropriate   probability 
paper    for   the   exponential   and   lognormal    distributions.       In 
some   cases,    plots    were  made    for    the    gamma    distribution. 

iii.      A   Chi -Square   Goodness-Of-Fit    test    was    per    - 
formed    for    lognormal    and   exponential    assumptions,    [Ref.    2] 

iv.      A   W    test,    due   to    Shapiro    and    tfilk,    [Ref.    3] 
was    used  to   test   the    lognormal   assumption    for    sample    sizes 
less    than   51,     (due    to    availability   of    the   tables), 
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v.      A   Kolmoqoro v-Smirnov    (K-S)    test    was    used    foi 
lognormal   and   gamma   assumptions    [Ref.    4]. 
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I I. STATISTICAL  CONSIDERATIONS 
A.   ROUGH  CHECK  METHODS 

1«   Boxplot  (Refs.  5  and  6] 

Purpose:  Many  batches  of  data  pile  up  in  the  middle. 
To  analyze  the  behavior  of  a  batch,  we  need  a  picture  of 
where  the  middle  (median)  lies  and  how  the  tails  relate  to 
it.   The  middle  is  generally  better  defined  than  the  tails 
so  we  want  to  see  more  at  the  tails.   Some  valuas  called 
outliers  are  so  low  or  so  high  that  they  seem  to  stand  apart 
from  the  rest  of  the  batch. 

One  technique  for  displaying  such  batchas  of  data 
which  is  often  mere  convenient  than  a  histogram  is  the  box- 
plot.   But  the  histogram  is  still  good  for  showing  the  shape 
of  the  distribution.   Boxplot  will  not  show  multiple  modes. 
There  may  be  several  "modas",  and  we  want  to  see  them. 

A  boxplot  is  obtained  by  first  calculating  the  lower 
and  upper  quartiles  and  the  median  of  the  batch  (sample)  of 
numbers  and  then  plotting  tha  numbers  en  a  horizontal  line. 

The  lower  quartiia  is  the  value  that  divides  the 
batch  into  two  parts,  with  1/4  of  ths  numbers  below  this 
value  and  3/4  above  it.   The  upper  qiartile  is  the  value 
with  3/4  of  the  observations  below  and  1/4  observations 
above  it. 
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The   next   step    is    to    draw   a    narrow   rectangular    bcx 
with    ends    corresponding   to    the   lower    and    upper    quartiles, 
and    to   display   the    median   point    by   a    plus    sign.      The    length 
of   this    box  is   called   the   i nterguart ile    range    or    H-Spread 
(Figure    1)  . 

To    the   above   are    added   lines    marking    1/2    H-Spread   on 
each    side   of   the   box.      Data    values   outside    the   lines     (out- 
liers)   are    marked    with   stars   and   rectangular   circles. 


***  1       +  i *      *  0030 

Interquartile  range 

1.5  interguartile  range 

2  interquartile  range 

-Figure  1  -  Boxplot- 

First,  the  data  are  ordered  such  that, 

X(1)  <  X(2) <  X(n)  where, 

X(1)     :    smallest    (minimum)    value 

X(n)     :    largest    (maximum)     valie 

media n=M    ,    center   of   batch 

if    n    is    an    odd  number,    H=X((n  +  1)/2) 

if    n   is    an   even   number    ,    M=  (X  (  (N/2)  +  1) +X  (  (n/2)  -  1)  ) /2 
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guartiles=H    ,    center  of    half   of    batch 
a  (H)=  (|  d  (M)  |  +  1)  /2  where       |  A|     means    integer   part   of    A 

inter guartile   distance  =    upper   quartile    -lower   guartile 

Here  are    good  boxplot   examples    for   the    exponential 
and   lognormal   model    (data   generated    randomly). 

BOXPLOT    EXPONENTIAL 

-I    ♦  I 0 


-Figure    2    -    Boxplot    Expor.er.tial- 
BOXPLOT    NORMAL    (log   is    taken    then    boxplot) 


-Figure    3   -   Boxplot   Normal- 

2 •      Hi s to  <jr a m   [ Re f .    5   and   6 ] 

Boxplot  is  good  for  comparison  but  a  histogram  is 
still  good  for  showing  shape.   Sometimes  histogram  is  more 
helpful  than  the  boxplot,  especially  for  detection  of  more 
bumps  in  the  data.   Histogram  prints  data  in  classes.   Exam- 
ples are  given  below  in  Figure  4. 
HISTOGRAM  EXPONENTIAL 

MIDDLE  OF    NUMBER  OF 
INTERVAL     OBSERVATIONS 


0. 

9 

*****  **** 

1. 

10 

*****  ***** 

2. 

2 

** 

3. 

4 

**** 

a. 

0 

5. 

0 

6. 

0 

7. 

0 

8. 

1 

* 

16 


HISTOGRAM  NORMAL  (log  is  taken  then  histogram) 

MIDDLE  OF    NUMBER  OF 
INTERVAL     OBSERVATIONS 


1.0 

1 

* 

1.5 

3 

*** 

2.0 

6 

*****  * 

2.5 

6 

*****  * 

3.0 

5 

***** 

3.5 

2 

** 

4.0 

1 

* 

4.5 

2 

** 

-Figure  4  -  Histograms- 

B.   STATISTICAL  DISTRIBUTIONS 

1  -   The  Exponential  Distribution 

The  probability  density  function  of  the  exponential 
distribution  is  [Ref.  2], 

f(x,A)=Ae     ,        x  >  0  and  X>  0 
Its  cumulative  distribution  function  is, 

F(x,X)=  YXe   wdr  =  i-e 

which   can   be   easily   evaluated. 

The   exponential  distribution    is    used   fr=guently   as   a 
time    to    failure    model    for    a    system    when    a   constant    failure 
rate    is   assumed   [Ref.    4],       But    for   r=pair    times,    it    can   be 
shown   that    it   is   not    an   appropriate    model.      Because   repair 
time    includes    diagnostic,    correction,    and    verification    tasks 
a  repair   time   density    must,    have    a    value    of    zero    at    time    t=0. 
It   should   then   increase   to    its    maximum    value   rapidly    and 
then    gradually    decrease   towards    zero    as    time   increases. 
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The  maximum  likelihood  estimator  of  tha  parameter  of 

the  exponential  distribution  from  given  data  is, 

n 
.L   Mi) 

A     1  1=1 

A  = where     x= 

Q 
X 

2.      Loqnormal    Distribution 

This   distribution   has   many    different   shapes   for   non- 
negative   variates.      It  is   skewed    to    the    right,    the    degree   of 
skewness   increasing   with   increasing   values   of   standard   devi- 
ation.     The   mean   and    standard   deviation    are    scale    and   shape 
parameters    respectively. 

A  random  variable  is  said  to  have  a  lognormal  dis  - 
tribution  if  the  logarithm  of  the  variable  is  normally  dis- 
tributed, 

0 

with   mean    and    variance 

2 

fJL   =   E    (lnX)         and  (J  =   Variance  (lnX) 

If    data    X(1),    X(2),     ,    X(r.)    and    if 

Y(1)=lnX(1),    Y(2)=lnX(2),    ,    Y(r.)=inX(n)    than, 

maximum    likelihood    estimators    of   the    lognormal    distribution 
are, 

A  A^    1  2 

JU=y  C= —     ,£,    (y(i)-y) 

n         i  =  1 
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th 
The   p        percentile    is  calculated   by  TRefs.    7    and    8] 

Median        Z  (.5)=0  X  =   e^ 

3    5 

th  JU +1.282  G" 

90    per.       Z(.9)=1.282  X  =    e 

0    9  r~ 

th  *  LC  +  1-6450 

95    per.       Z(.95)  =  1.645  X  =    e^ 

D.95 


3 .   Gamma  Distribution 

k   random  variable  x  has  a  ganma  distribution,  and  it 
is  referred  to  as  a  gamma  random  variable,  if  and  only  if 
its  probability  density  is  given  by,  [Ref.  9] 

|    1       a-1   -x/b 

j x     5  for  x  >  0 

f(x)=  I   a 

b    Ha) 

0  elsewhere 

where     a  >  0  ani  b  >  0 

a  is  shape  parameter 

b  is  scale  parameter 

"Pis  the  gamma  function 

The  mean  and  variance  of  the  gamma  distribution  are 

given  by, 

2 
mean=a.b  and  variance=a.b 

Special   cases    of    the    gamma    distribution    play    an    im- 
portant   roles   in    statistics.      For   instance,    for    a=1    and   b=Q, 
we   get   the    exponential  distribution    with     A=1/Q. 
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For  the  gamma  distribution  the  scale  and  shape  para- 
meters can  be  estimated  from  data  in  many  ways;  here  are  two 
a.   Method  Of  Moment  (M.O.fl)  estimators 


2 

(sample   mean)  sample    variance 

a= b  = 

sample    variance  sample    mean 

2 
where  mean  =  a.b        and   vanance=a.  b 


b.       Maximum   Likelihood   Method"     (M.L.E.)    estimators 
This   method   is    more    complex   as    compared   to    the 
M.O.M.;    it    requires    numerical   calculation.      There    is    an    APL 
function    named    ESTGAM    that    gives   directly    M.L.S.     (App.C). 
It    is    used    because    it    provides    a    batter    estimator    of    the   pa- 
rameters  a    and   b    in    large   samples,    provided   the   data    truely 
come    frcm   a   gamma    model. 

C.       PROBABILITY    PLOTTING 

Plotting    data    points    on    probability    paper    is    quite    easy 
and    does  not   require    many   calculations.       It    is   a    visual   exa- 
mination  rather    than    statistical   calculation    (Figure    5) . 

What   can    we    learn    from    a    probability   plot? 

1.  Visual   indication   of   assumed    iistribu-ion . 

2.  Examine    departures    from    assumed    distribution. 
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Jty 


J 


-1 


figure    oa 


ir'igara    bD 


5a.   Tails  longer  than  normal 
5b.   Tails  shorter  thin  normal 

-Figure  5- 


1 .      Lognormal   Distribution 

Consider   the    order    statistics    of   a    sample 

x(1)     <   x(2)     <    <    x(N)     Then,    x  (i) 

is  an  estimate  of  the  (i/N)  th  percentile  of  X 
If  X  is  N(0,1)  then  we  expec4: 


-1 
x(i)^r^      (i/N) 


for        1*1,2,3 

-1 


Then   a    plot   of  X  (i)     vs  jS      (i/N)     should   give   a 
straigh-    line.      The   only   problem    is 


-1 
$      (i/N) -infinity,    so  we   us=     j8T   (i/(N+s)),    s»1. 

"1E7- 


-1 

r 

PROCEDURES 


a.  Use   normal    probability    paper 
Plot    (i/N+s)    vs    X(i) 

b.  Use   norm plot   in    APL   workspace   [Ref.    21] 
Use    UINITAB 

NSCORES    of    Cl,put    in    C2 

PLOT   C2    vs    C1 


NSCORES    are     (i/N  +  s) 
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Consider   log  of    data   Y(1)    <   Y  (2)     < 

<   Y(N)       where    Y(i)=LnX(i)     and    then    apply    above 

procedure   no    get    lognormal    plot. 
2.      Exponential   Distribution 

Suppose    X    is    exponential    A     with   sample    size    ?. ,    so 
Pr(X    <   x)  =  1-exp  (-Ax)  where   X    >    0.       Then    (i/r.  +  s)    is 

an   estimate   of   Pr  (X   <    X  (i)  )  =1-exp  (-Ax  (i)  ) 

So    if   the    data   really  oome   from    an   expoiential   dis- 
tribution   (i/N+s)       vs.      1-e  xp  (-Ax  (i)  )  ,    should   yield   a 
straight    line. 

PROCEDURES 

a.  In  (i/N+s)    vs.    X (i)     oi    rectangular    paper, 

b.  EXPONPLOT    in    APL    workspace,    [Ref.     21] 

c.  MINITAB, 

d.  Use   Gamma   probability    paper    for    ALPHA=1 
BETA=Q   where   Q   is    exponential   parameter    so    that, 

f  (x)=  (1/Q)  *exp(-X/3)     for    X   >    0 
3  .      Gamma   Distribution 

There    is   an   APL    program   function    named   3AM    that    uses 
estimated   gamma    parameters    and   gives    test    statistic    D.      The 
program    also    gives    a    matrix    for    probability    plotting    pur- 
poses,      (see   Appendix    C) 
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D.   STATISTICAL  TESTS  FOR  DISTRIBUTIONAL  ASSUMPTIONS 

We  need  to  know  whether  we  are  actually  dealing  with  a 
sample  from  an  exponential  population  or  whether  the  data 
values  of  random  variables  have  a  lognormal  or  some  other 
distribution.   A  statistical  test  of  significance  of  a  dis- 
tributional assumption  provides  an  objective  technique  for 
tentatively  assessing  whether  or  not  an  assumed  model  pro- 
vides an  adequate  description  of  the  observed  data. 

There  are  three  basic  steps  involved  in  statistical  test 
methods. 

1.  A  test  statistic  is  calculated  from  the  data, 

2.  The  probability  of  obtaining  the  calculated  test 
statistic  is  determined, 

3.  Assessment  is  mads  of  The  adequacy  of  the  assumed 
distribution. 

a.  If  the  probability  of  obtaining  the  calculated 
test  statistic  is  "LOW",  one  can  say  that  the  assumed  dis  - 
tribution  dees  not  provide  an  adequate  representation. 

b.  If  the  probability  of  obtaining  the  calculated 
test  statistic   "HIGH"  ,  then  the  data  provides  no  evidence 
that  the  assumed  distribution  is  not  adequate,  judged  by  the 
test. 
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The  definition  of    "LOW"    or    "HIGH"    depends    on    the 
user's   preferences,    and   the    consequences    of   rejecting    the 
distribution.      For    example,    ths    rejection    statistic   is   typi- 
cally   be   0.01,    0.05,    or    0.1.      For   the    purpose    of    this    study 
0.05    was   selected    as    the   reject    criterion. 

It    should    be    pointed   out    that    a   statistical   test, 
although   it   allows    one  to  reject    an    assumption   as    inadequate 
does   not   prove   that    the   assumed    distribution   is   correct. 
1 •      The  Chi-Sguared   Soodness-of ;Fit    Test 

This   test    is    one    of    the    oldest    and    most    commonly 
used    for   evaluating   distributional   assumptions.       Basically, 
the    given    data    are    grouped    into    cells    and   compared    to    the 
expected   number    of    observation    in    the    cells    based    on    the 
assumed    distribution.      Then    if    the    test    statistic,    calcula- 
ted   from  this   comparison    exceeds    a   critical   chi-square    value 
the    assumed   distribution    is    rejected. 

The  only    problem    is    dividing    the    data    into   cells. 
If   the   number   of   observations  is   small    it    is   usually    sugges- 
ted   to   use    a    number   of  cells   as    large    as    possible,    subject 
to    the   restriction    that    it    should   not    exceed   N/5     (   N    is   sam- 
ple   size). 


2 '4 


The   computations   involved   in    the    Chi-Squared    Good   - 
ness-Of-Fit    test    were    made    by  a    compiter    program,    in    which 
are    given   as    inputs   the   assumed   cumulative    distribution 
function,    the   number    of   observations,    the   number    of   equipro- 
babie    cells,    and   the    number    of    parameters   estimated    from   the 
sample.      The   outputs    are   the   Chi-Square    statistic    and    its 
probability   of   exceeding   a    chi- square    value   for   a    given   num- 
ber   of   degrees    of    freedom. 

The   procedure    used    is  described    as    follows:^ Ref .    2] 
a.      The    cell    boundaries    are   determined   from   the 
assumed    cumulative    distribution    as    the    values   such   that    the 
estimated    probability    of    the   observation    value    falling 
within   a   given    class    is    1/K .      For    saoh    class, 

Pr  (X<x  (i)  )    =    i/k,    where 

x=the  random  observation  to  be  assigned  to  xhe 
i-th  cell 

x  (i) =the  i-th  sell  boundary  to  be  solved  from 
the  above  formula 

k  =  the  number  of  cells,  in  this  study  k  =  n/5 

The  lower  bour.d  of  the  first  c  =  ll  and  -he  upper 
bound  of  the  last  cell  are  the  smallest  and  the  largest 
values  cf  the  observations. 
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b.  The  expected  number  of  observations  for  each 
cell  E(i)  is  equal  to  n/k  for  each  sell. 

c.  The  number  of  observed  values  in  each  cell  M(i) 
is  counted  based  on  the  results  of  the  above  equation. 

d.  The  test  statistic  is 


2  k   t    K         .    2   T 

n     |    1=1  I 


2 

e.      The   computed    value      ^    '  is    used   to   compute    the 

level   of  significance,    or   the    probability   of   a   3hi-Square 

value    with    deqrees    of    freedom   equal   to    K-M-1    where    H    is   the 

number  of   estimated    parameters   from    the    sample   to    exceed 

calculated    chi-square 

2  2 

<X=Pr(         (V)    >    X   )  with  V   degrees   of    freedom 

1-k 

If  oC   is    less    than   or  equal    to    .05    the    assumed    dis- 
tribution  can    be   rejected   as    inadequate. 

2-      The    Kolmogoro  v-Smimov    Test    (K-S    Test) 

K-S   test    examines   the   cumulative    distribution    func- 
tion   rather   than   probability   density    function    .       It    simply 
compares  the   sample    cumulative   distribution    function    S(.) 
to   the   hypothesized   distribution   function    ?(.).       As   a    mea- 
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sure    of  comparison    the  test    uses,    D    =    Max | F  (. ) -S ( .)  |  ;    this 
is   the   main    difference   between    the   Cai-Squared   3.0.?.      test 
and    the    K-S.    test.      [ Ref .    U  ] 

On    the   other   hand   here   we   ara    dealing   with   Type    I 
error,    that   is   the   hypothesis  is   true    and  reject    it.      The 
probability  of   commit inq    a    Type   I    error    is    called    the    signi- 
ficance   level   for    the    test    alpha.      Here    alpha    is    taken    to    b- 
0.05.      This    means    no    mere   than   a   5    %    chance   of   rejecting   th- 
true    hyphothesis. 

X    has    CDF    F(x)=    Pr  ( X    <    x)  , 

Take   a   sample    X(1)     <   X(2)    < <   X(N), 

We    want    to    estimate        F(x)=    Pr  (X    <    x)  =   S     (x) 

n 


S     (x) 

n 


■    0 


1/N 
2/N 


K/N 


for    X  <    X  (1) 


"       X    >    X  (1)  .  AND.X    <    X(2) 
"       X    >    X(2)  .AND.X    <    X(3) 


"       X   >    X(K)  .  AND.X    <    x  (K+1) 


1  X   >    X (N) 

Find   D=    MAX  |     S(x)  -F  (x)  | 

Using   Table   9    of    reference   4    for    given   eignificar.ee 


levex, 


If  D  exceeds  the  tabulated  value  that  comes  from  the 


table,  reject  the  hypothetical  distribution. 
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If    D   does   not    sxcz^Z   the   tabulated    value    that    comas 
from   table,    do   not   reject   the  hypothetical    distribution. 

There   are   two    kinds    of   tables    for   the   K-S    test.      One 
of   these   is   used    if   the   parameters   are    estimated    from    data, 
and    the   other   one    is    used   if   the    parameters   are    given.       Here 
parameters    are    estimated    from   data. 

3.      The    W   test    for  Lo^normal    £ssamp_tion 

The  1  test  is  shewn  in  Ref.  3  - o  be  an  effective 
procedure  for  evaluating  the  assumption  of  normality  against 
nonnormal  alternatives,  even  if  only  a  relatively  small  num- 
ber of  observations  are  available.  Hahr  and  Shapiro  suggest 
that  the  W  test  may  also  be  used  to  evaluate  tha  assumptions 
of  the  logncrmal  distribution.  This  follows  beoause  of  ^he 
property  that  if  the  logarithms  of  the  observations  follow  a 
normal  distribution,  then  the  original  values  of  the  obser  - 
vations   are    lognorraally   distributed. 

The    following    procedure    was    applied   usi?.g    tables 
from    Reference    2. 

a.      The    observations   are   ordered. 


b. 


X(1)    <   X(2)     <   ...<   X  (i) <    X(n-1)     <   X  (n) 

The    following   parameters    are   computed, 


S    =    ^[lnx(i)]      -    -i 


X     lnx(i) 
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ii.      If   n    is    5ven,    k=n/2 

n    is    odd,    k-  (n-1)  /2 

B=    .g-C    a(k-i+1)«C  (lnx(n-i+1))-(lnx(i))  ]] 

where  tha    values   of    a(k-i+1)    for    i=1,2,..k 

are    given    in    table   IX   [Ref.     2]   for   n  =  3    to   n  =  50  • 

iii.      The    test    statistic    W    is 

2        2 
W=    B  /   S 

c.   The  approximate  probability  of  obtaining  the 
calculated  value  of  W  can  be  obtained  from  the  formula 

v         i w-e  i 

Z=A+77ln  J— -| 
,f         h-wj 

using  the  values  of  A  ,77  and  6  given  in  Table  XI  of  Ref.  2 
for  the  appropriate  sampla  size  and  then  using  standardized 
normal  distribution  to  determine  the  probability  of  obtain- 
ing a  value  less  than  or  equal  to  Z,  which  is  the  signifi- 
cance level  of  the  test. 

q  =  Pr(Z   <  Z) 

q 

If  q  is  less  than  or  equal  to  .05/  the  selected  le- 
vel of  significance  in  this  study,  the  logr.ormal  is  rejected 
as  an  inadequate  assumption. 
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III. DATA    ANALYSIS    RESULTS 
A.       RESULT    OF    TESTS    FOR    DISTRIBUTIONAL    ASSUMPTIONS 
1 .      Summary   of   results 

Tables    2,3    and  4    summarize   the    results   of    the    sta   - 
tistical   +est   analysis.      In    all    cases    the   5    %   level   was   cho- 
sen   as  the   reject    criterion.      The   46    sets   cf   data    used   in 
this    study    can    be   divided   into  three    parts    as   indicated   in 
Section   on    IB2. 

The   following    is    an    example    of    the   computer   program 
summary    table   which   includes   the    results    for    the    Chi-Square 
test    and   calculated    parameters   from    the    sample   data. 
SET    NO    18    -    AN/GXS-2  (V) 

SAMPLE    SIZE      N    =       26  NO.     OF    CELLS      K    =       5  (a) 

SAMPLE    MEAN    =       19.70  STANDARD    DEV    =       22. b2        (b) 

EXPONENTIAL  LOGNORMAL  DIFFERENCE 


PARAM1 

0. 

05 

2. 

,56 

(c) 

PARAM2 

0. 

,78 

MTTR 

19. 

70 

19. 

,  13 

2. 

.96 

1 

(3) 

50- 

TH    PERCNT 

13. 

65 

12. 

,93 

5. 

.61 

% 

90- 

TH    PERCNT 

45. 

35 

40. 

,22 

12. 

.77 

% 

95- 

TH    PERCNT 

59. 

00 

55. 

,46 

6. 

,39 

% 

CHI 

-SQR    STAT 

5. 

92 

0. 

.54 

<e) 

DEG 

OF    FREED 

3 

I 

; 

) 

(f) 

SIGNIF    LEVEL 

0.  115E+ 

00 

3.764E+00 

(g) 

30 


r  able  2 

'OTAL  DATA  ANALYSIS  RESULTS 


S=t  Sample 
No    Size 


1 

43 

2 

43 

3 

30 

i» 

43 

5 

33 

6 

21 

7 

59 

8 

20 

9 

49 

10 

50 

11 

50 

12 

43 

13 

30 

14 

26 

15a 

33 

16 

30 

17 

50 

18 

26 

19 

57 

20 

57 

21 

157 

22 

28 

23 

36 

24 

58 

25 

20 

26 

92 

27 

21 

28 

33 

29 

47 

30 

35 

31 

81 

32 

22 

33 

45 

34 

171 

35 

31 

36 

37 

37 

60 

38 

20 

39 

95 

40 

22 

41 

35 

42 

51 

43 

38 

44 

87 

45 

22 

46 

47 

E    X 

P    . 

LOG    NORMAL 

GAMMA 

TE 

ST 

TEST 

W    K 

-S 

TEST 

Pr. 

Ch. 

3o. 

Pr. 

Ch. 

M.     M. 

PI. 

Sq. 

Rej. 

PI. 

PI. 

Sq. 

rp 

T. 

Rej. 
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B 

R 
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G 
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A 

R 
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B 
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+ 
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B 

R 

R 

R 

+ 
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B 

R 

+ 

B 

? 
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R 
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* 
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A 

B 

B 
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R 

p 
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3 
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♦ 
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F 
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— 
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G 

R 
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P 

F 

R 

R 

R 
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P 

R 

•f 

G 

j 

A 

A 

B 
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P 

R 

+ 

F 

F 

R 

A 

a 

* 

R      R 

P 

R 

+ 

G 

G 

A 

A 

A 

A       A 

3 

R 

♦ 

F 

F 

R 

A 

R 

* 

A      A 

B 

R 

+ 

P 

p 

R 

R 

R 

♦ 

A      A 

P 

R 

+ 

G 

j 

A 

A 

A 

A      A 

G 

R 

+ 

P 

n 

A 

A 

R 

R      R 

B 

R 

+ 

B 

s 

R 

R 

R 

+ 

R      R 

P 

R 

+ 

G 

J 

A 
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A 

R       A 

F 

A 

F 

F 

A 

A 

A 
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B 

R 

♦ 

B 
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A 
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+ 
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B 

R 
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F 
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R 
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P 

R 

+ 

G 

F 

R 

R 

R      R 

F 

A 

G 

j 

A 

"I 

& 

R      A 

B 

R 

♦ 

P 

B 

R 

R 

A 
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R 
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B 
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R 

R 

+ 
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F 

R 
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G 

B 

R 

'J 

A 
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F 

R 

+ 

F 
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R 
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B 

R 
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B 

B 

R 

"I 

R 

* 
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B 

R 

+ 

G 

j 
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A 

A 
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P 

R 
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G 
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A 

A 

A 
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3 

R 

+ 

F 

j 
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A 
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P 

R 

•f 
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R 
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R 
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R 
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R       R 
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B 
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+ 
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A 
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R       R 
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+ 
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F 
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R 

* 

A       A 

F 

R 
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F 
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+ 

G 

G 

A 

~£ 

R 

A       A 

B 

R 

♦ 

B 

3 

R 

R 

R 

♦ 

R       R 

2 

Criteria    for    rejection   Pr(X    )     cr    Pr  (W)    less   than    .05 

or    D   value    for  K-S    test    bigg=r    than   tabulated    value. 
*      To    be  determined    from   probability   plot    or    BoxDlot 
where  B=3ad,    P=Poorr    F=Fair,    G=Good 

A=Do   not    reject,    R=Reiect 
Pr.Pl.=    Probability   Plot,  Bo.Pl.=    Boxpiot 
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LOGNORHAL    AND   E 

pifl 

["Sample- 

-I — b-x  '  P  'J  " 

SO     1 

1 

Size 

1       P(X    ) 
j 

1- 1 

1 — «ra — 

ur*J7T*g~ 

2 

43 

0.00166 

3 

30 

0.48E-4 

4 

43 

0.59E-7 

5 

33 

0.018 

6 

21 

0.245 

7 

59 

0.2  38E-6 

8 

20 

0.0101 

9 

49 

0. 119E-6 

10 

50 

0.5  96E-7 

11 

50 

0.00011 

12 

43 

0.596E-6 

13 

30 

0.28E-5 

14 

26 

0.  129E-4 

15 

50 

0.208E-4 

15a 

33 

0.216E-4 

16 

30 

0.8  11 E-4 

17 

50 

0.0323 

18 

26 

0.  115 

19 

57 

0.596E-7 

20 

57 

0.0 

21 

157 

0.0 

22 

28 

0.381 

23 

36 

0.808E-4 

24 

58 

0.  132E-3 

25 

20 

0.0334 

26 

92 

0.6  56E-6 

27 

21 

0.  174E-2 

28 

33 

0.018 

29 

47 

0.69E-2 

30 

35 

0.254 E-4 

31 

81 

0.  186  E-  3 

32 

22 

0.468E-2 

33 

45 

0.427E-3 

34 

171 

0.5  96E-7 

35 

31 

0.456 

36 

37 

0.0258 

37 

60 

0.542E-5 

38 

20 

0.0334 

39 

95 

0.298E-6 

40 

22 

0.0716 

41 

35 

0.0476 

42 

51 

0.193 

43 

38 

0.470E-2 

44 

87 

0.182E-2 

45 

22 

0.468E-2 

46 

47 

0.7  57S-5 

Table    3 
•ONENTIAL    RESULTS 


2                                                K  —  S  ^^s "** 

P(X   )          P(tf)          |  Ts.St.  Tab.^Va.l 

"TJ7B39" ' 073*35 ' 7T4T 7115"" ' 

0.358  0.025  .079  .135 

0.0588  0.0745  .205  .161 

0.0005  0.0  .238  .135 

0.0128  0.0328  .148  .154 

0.0468  0.0002  .27  .193 

0.596E-7  .172  .115 

0.0143  0.017  .235  .190 

0.459  0.31  .217  .126 

0.0445  0.099  .212  .125 

0.979  0.632  .086  .125 

0.0421  0.19  .143  .135 

0.0186  0.01  .3  .161 

0.292  0.238  .151  .174 

0.00408  0.001  .245  .125 

0.484  0.06  .211  .154 

0.0321  0.0  .199  .161 

0.494  0.583  .125  .125 

0.764  0.36  .084  .174 

0.360  .23  .117 

0.116E-4  .441  .117 

0.0  .073  .07 

0.526  0.624  .145  .167 

0.173E-3  0.0074  .132  .147 

0.0255  .143  .116 

0.912E-2  0.35  .16  .190 

0.380E-4  .102  .092 

0.231E-2  0.175  .301  .193 

0.554  0.345  .143  .154 

0.711  0.524  .106  .130 

0.592  0.617  .263  .149 

0.861E-4  .121  .098 

0.063  0.59  .219  .189 

0.315E-3  0.013  .147  .132 

0.596E-7  .069  .068 

0.724  0.97  .117  .159 

0.571  0.11  .083  .145 

0.430E-2  .123  .114 

0.094  0.896  .099  .110 

0.22  .091  .091 

0.201  0.424  .156  .189 

0.126  0.48  .125  .149 

0.103  .092  .124 

0.242  0.04  .265  .144 

0.193  .106  .095 

0.340  0.57  .224  .189 

0.192E-3  0.004  .179  .13 
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Table    4 

GAMMA    DISTRIBUTION    TEST    RESULTS 
(Kolmogorov    Smirnov  G.O.F.Tast) 


Set 

D   value   of 

A/R 

D   value   of 

A/R 

Tabulated 

No 

M.O.M.Est. 

_  — — 

M.L.E.Es-. 

Value 

1 

.36  2 

R 

.  154 

R 

.135 

2 

.286 

R 

.  182 

R 

.  135 

3 

.22 

R 

.  185 

R 

.161 

4 

.45  8 

R 

.  344 

R 

.  135 

5 

.37 

R 

.  21 

R 

.154 

6 

.23 

R 

.  164 

A 

.  193 

7 

.29  9 

R 

.  138 

R 

.115 

3 

.  18 

A 

.  154 

A 

.190 

9 

.09  3 

A 

.  11 

A 

.126 

10 

.  143 

R 

.  147 

R 

.  125 

11 

.08  1 

A 

.081 

A 

.125 

12 

.124 

A 

.  107 

A 

.  135 

13 

.  148 

A 

.  152 

A 

.161 

14 

.  168 

A 

.  153 

A 

.  174 

15 

.  158 

R 

.  162 

R 

.125 

15a 

.17 

R 

.  114 

R 

.  154 

16 

.  23  1 

R 

.  162 

R 

.161 

17 

.  136 

R 

.097 

A 

.125 

18 

.23  2 

R 

.  143 

A 

.174 

19 

.07  3 

A 

.074 

A 

.117 

20 

.  145 

R 

.  109 

A 

.117 

21 

.192 

R 

.  163 

R 

.07 

22 

.12 

A 

.  112 

A 

.157 

23 

.285 

R 

.  20  1 

R 

.  147 

24 

.26 

R 

.  164 

R 

.116 

25 

.  172 

A 

.  155 

A 

.190 

26 

.22  2 

R 

.  14 

R 

.092 

27 

.  186 

A 

.  187 

A 

.  193 

28 

.  172 

R 

.09 

A 

.154 

29 

.  187 

R 

.113 

A 

.130 

30 

.  1 

A 

.098 

A 

.149 

31 

.086 

A 

.085 

A 

.098 

32 

.  105 

A 

.114 

A 

.189 

33 

.  185 

R 

.  142 

R 

.132 

34 

.284 

R 

.  13 

R 

.063 

35 

.  153 

A 

.098 

A 

.  159 

36 

.268 

R 

.  153 

R 

.145 

37 

.286 

R 

.  172 

R 

.  1  14 

38 

.  124 

A 

.  107 

A 

.19 

39 

.336 

R 

.  161 

R 

.091 

40 

.08 

A 

.085 

A 

.184 

41 

.  154 

R 

.  099 

A 

.  149 

42 

.  163 

R 

.099 

A 

.124 

43 

.  13  1 

A 

.  126 

A 

.  144 

44 

.08  8 

A 

.  075 

A 

.095 

45 

.  12 

A 

.08  9 

A 

.  189 

46 

.23  3 

R 

.  175 

R 

.13 
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Notes : 

a.  The  number  of  equiprobabLa  cells,  X,  was  chosen 
as  N/5  where  N  the  sample  size. 

b.  The  sample  mean  and  standard  deviation  are  cal- 
culated based  on  the  maximum  likelihood  estimates. 

c.  PARAM1  for  the  exponential  distribution  is  the 
reciprocal  of  the  sample  aean,  for  the  lognormal  distribu- 
tion  PARAM1  and  PARAK2  are  jj  and  G ,  the  parameters  of  xhe 
lognormal  distribution. 

d.  The  MTTR  and  the  50th,  93th  and  95th  percentiles 
are  calculated  from  the  sample  and  are  based  on  the  rela- 
tionships between  the  calculated  parameters  and  their  dis- 
tribution functions.   The  percentage  difference  is  between 
the  exponential  and  lognormal  MTTR  and  percentiles. 

e.  The    Chi-Square    statistic. 

f.  The   number   of   degrees   of    freedom    (K-M-1)     is    K-2 
for    the    exponential    distribution   and    K-3    for   the    lognormal. 
This    is    because    one    parameter   is   estimated    from    the    sample 
in   the   exponential    case   and    two    in    the    lognormal    case. 

g.  The    level    of    significance    is   the    probability    of 
a  Chi-Square   variate    with   the   specified    degrees    of    freedom 
exceeding  the    calculated   Chi-square   statistic. 


34 


For   the    exponential    model,    only    the    Chi-Squared   test 
was    used.       In   6    oat   of  46   cases    we   would   accep-    the    exponen- 
tial   model.      If    W3   compare    the   exponential    vs.    lognormal    mo- 
del   using   the   Chi-Sguare    test,    test    prefers    the    lognormal 
model   to   the   exponential    model.       In   only   one   case    would   the 
exponential   model   be   accepted  and   the    lognormal   rejected. 

exponent ial 
A  B 

A  I     5       I     1  9    |  2'4 

lognormal  j j | 

R  |     1       |     2  1     |  22 

6  40 

Therefore    the    lognormal    model    is    a    better    descriptor 
with    respect    to   the    exponential    model. 

In    most   cases   the   lognormal   assumption   works    reason- 
ably   well.      In   those    cases    in   which   the    results    of   the   chi- 
sguared    test,    W    test    and    K-S    test    do    not    agree,    a    plot    on 
lognormal   probability    paper    was    used    to    determine    the    appro- 
priateness   of   this    assumption    (*    in    Table    2).      Comparison 
between    test   methods    were    lade   in   three    ways    for    the    lognor- 
mal   model   and   tabled    below. 

a.       Chi-Sg.    vs.    K-S    test     (46    cases) 

K-S 

A  R 

A         |     15    |    10    J         25 
Chi-Sg.  ,----!----, 

R  |       3    j     1  8    |  21 

18         23 
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b.  Chi-Sq.    vs.    W   test    (34   cases) 

W 
A  R 

A  |     19    |        2    |  21 

Chi-Sq.  J | | 

B  |       4     |        9     |  13 

23         1  1 

c.  W   vs.    K-S    test.     (34   cases) 

K-S 

A  B 

A         I     13    I     10    |  23 

B  j       3    |        8    |  11 

16         1  8 

For   the    lognormal  methods   the    test   methods   sometimes 

give    different   results.      The    K-S    test    and   the   chi-square 
test    agree    in    74    %   of    the   cases.      Th=    chi-square    agrees    with 
the    9   tesr    in   82    %    of    the   cases.       A    and    K-S   test    agree    62    *. 
Therefore,    there   is   no  significant    difference   among   the 
three   test    methods. 

A    previous    study    has   been   done    using    24    data    sets    of 
electronic    systems    and  equipments,    using    only    the    H    and   chi- 
squared   tests    for   the    lognormal    model    [Ref.    1].      The    table 
below    gives   a   comparison    between   the    exponential    and    lognor- 
mal   for   the   previous    study.       For    the    exponential    model   do 
not    reject     (A)    or    reject    (R)    decision    was    made   based    on   the 
chi-squared   test.      For   the    lognormal    model    accept    by    either 
the    chi-square   or   W    test    was    used   as    a    criterion. 
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lognor  mal 


a 


ex  ponential 
A           B 

1     7 

1 

1     1 

14    | 

"31 

21 
3 

8 

1  6 

The  xables  below  show  same  kind  of  results  for  the 
current  study.   Only  difference  between  the  two  studies  is 
there  are  different  kind  of  systems  and  equipments  (e.g. 
mechanic  or  electronic)  in  the  current  study. 


ex  ponential 
A            R 

logncr mal 

A 

R 

I    5 

23    | 

"31 

23 
18 

6 

ao 

For  all  cases  (46  cases) 

exponent  ial 
R 


lognormal 


I  1   I   7  |     3 

1     1  2 

For  electronic  eqaipments  only  (13  cases) 


These  results  indicate  that  the  loanormal  model  is 
better  than  the  exponential  model  for  repair  time  data. 

For  the  gamma  model,  the  MLE  estimator  is  more  po  • 
werful  than  the  MOM  estimator,  with  the  K-S  test  used  for 
testing  the  data.   The  decision  was  made  to  accept  the  re  ■ 
suits  of  the  ML2  when  the  two  estimators  gave  different 
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results,  since  in  no  case  did  the  MOM  suggest  acceptance 
when  the  MLE  said  to  reject.   In  the  eight  cases  where  there 
was  a  difference,  therefore,  the  MLE  was  used  as  the  crite  - 
rion.   This  gives  the  following  results  for  the  gamma  model. 

MOM 

A     R 

A    |  19  |   8  |    27 

MLE       J | | 

R       0  |  1 9  |    19 

19    27 

Therefore  using  the  MLS,  we  would  accept  the  gamma 

model  in  27  out  of  46  cases.   In  16  of  these  27  cases,  the 
lognormal  model  is  also  accepted  and  in  5  additional  cases 
all  three  distribution  assumptions  would  be  accepted.   Again 
only  one  case  (case  6)  would  we  accept  the  gamma  and  expo- 
nential model  but  reject  the  lognormal. 

We  can  also  compare  lognormal  and  gamma  model  using 
the  K-S  test. 

gamma 
A     R 

A    I  13  I   5  I    18 

lognormal       I j 

R    |  13  I  15  |    28 

26    2  0 

The  above  table  shows  that  the  K-S  test  prefers  the 

gamma  model.   27  out  of  46  cases  accepted  the  gamma  model, 

but  only  18  out  of  46  cases  accepted  the  lognormal  model. 
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Tables  6-9,  (Appendix  B)  sumaariza  the  results  for 
all  cases  and  for  each  of  the  three  different  types  of  data. 
Table  10  (Appendix  B)  gives  total  results  in  terms  of  sample 
sizes  of  20-29,  30-39,  40-50  and  over  50.   For  exponential 
assumption,  increasing  sample  size  increases  rejection  per  - 
centages.   For  gamma  assumption,  increasing  sample  size  also 
increases  rejection.   In  lognormal,  increasing  the  sample 
size  decreases  rejection  up  to  sample  sizes  of  forthy,  then 
rejection  increases  rapidly. 

Because  of  the  large  number  Df  data  sets,  only  some 
interesting  cases  of  the  data  analysis  and  probability  plots 
for  those  sets  are  discussed  in  the  next  section. 
2-   Discussion  of  speci  fie  cases 

The  following  discussions  illustrate  some  cases 
which  are  of  special  interest.   The  data  analysis  results 
and  plots  are  given  in  Appendix  A. 

a.   Set  No.  18  -  AN/GXS-2(7) 

This  case    is   an    illustration    where   ail   three 
models    are    accepted    (not    rejected)    for    all    tests.       Data    were 
collected   from    26    maintenance   tasks.      In    this   case,    as    in 
most    other    of   the   formal   maintainability    demonstration    (ca- 
ses   8-20) ,    -he  tests    were    performed    in    accordance    with   the 
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U.S.  military  standard  for  maintainability  demonstration 
(MIL-STD-471)  [Ref.  23]  in  which  sample  selection  was  ge  - 
nerated  from  knowledge  or  estimates  of  component  failure 
rates  and  tests  were  conducted  using  service  technicians  of 
the  appropriate  skill  levels  and  given  prior  training  on  the 
equipments. 

This  case  fits  the  lognormal,  exponential  and 
also  gamma  for  both  estimators.   All  three  tests  give  high 
significance  for  lognormali ty. 

b.  Set  No.  6  -  Reactor  Coolant  Pump 

This  case  is  the  only  one  the  exponential  model 
accepted  and  the  lognormal  model  rejected.   However  the 
plots  are  not  good.   Further  investigation  is  necessary  in 
this  case  but  since  only  the  data  was  available,  it  was  not 
possible  for  us  to  do  so. 

c.  Sets  No.  28  and  41  -  Main  Rotor  Head  and  Blades 
These  two  cases  indicate  a  very  good  example  of 

logncrmality.   All  three  tests  signal  "do  not  reject"  fcr 
the  lognormal  assumption.   The  plots  are  also  very  good. 
The  gamma  family  does  not  well  represent  the  data  sets. 

d.  Sets  No.  1  and  2  -  French  Fessenheira  Pumps 
(  Down  /  Repair  )  times 
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These   two    cases    are    discussed    as   unsuccessful 
for    gamma  and   exponential   but    good   for    lognormal    since    ax 
least    2    out  of   3    tests  give    acceptance    for   the   lognormal 
assumption.      For   repair   times   the   W    test    rejects    the    assump- 
tion   while    the  Chi-squared    and   K-S    tests    do   not,    but    down 
times    are   rejected    by    only    the   K-S    test.      There   are   some 
outliers   around    5%,    but    the    rest    of    the    data   shDw    a    good 
lognormal    fit    in   the    probability   plot. 

The    existence   of   outliers    could   be   a    reason   for 
rejection   by    one   of   the    3   tests.      For    downtime   data,    after 
trimming   the   two   high   outliers   the   lognormal   assumption   was 
also    accepted   by   the    K-S    test.      For    repair    time,    after   trim- 
ming   one   outlier   the   data   shows    high    significance    (0.15)     for 
the    W    test.      Thus,    in    some    cases   it    would   be   valuable   to 
examine    in    more   detail  the    reason    for    the    high   or    low   outli- 
ers.     In    this   case,    only    the   data   was    made    available    so    that 
this    was   not    possible. 

e.      Set    No.    4    -    Heater   Drain    Pump 

Set    No.    5    -    Main   Feedwater    Pump 

Set  No.  4  illustrates  the  case  where  all  three 
tests  for  the  lognormal  model  are  rejected.  Here  outliers 
are    important    factors.      The    effect    of    outliers    can    be    seen 
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easily   from   the   boxplot,    histogram   and    probability    plot.      No 
5   is    almost   the   same    as   no    4,    the   only    difference    between 
the    two   being   that    the  K-S    test    barely    accepts   the    lognor   - 
mality   in    set   no   5.       For    these   cases,    trimming   the    outliers 
did    not   change   anything   in    so   far    as    significance    levels    or 
acceptance   or   rejection  are    concerned. 

f.      Set    No.    9    -    AN/ARC-162    (Intermediate) 

Set    No.    10   -    AN/ARC-163    (Intermediate) 

These    two   cas^s    are    discussed   together,    because 
they    give   opposite    results.       Derailed    reports   exist    in    both 
cases.      Test    technicians    were   used    for    running   the    tests    and 
removina   the    faults.      [Refs.    12    and    13] 

The   exponential    assumption    is    rejected   in    bo*h 
cases.    In   set   No   9    the   lognormal    assumption    is   rejected   only 
by   the   K-S    test    while    the   other    tests    give    high   signifi- 
cance.     The   probability    plot   and    boxplot    also   signal    that 
"do    not   reject"    is    appropriate    for   this   case.      The   gamma 
model    fits    this    set    for   both   estimators. 

In    Set    No.    10    only    the   'A    test    accepts,    but    the 
Chi-Sguared   test    barely    rejects    (0.0'45)     and  the    K-S    test 
rejects    the   lognormal    model.      However    the    boxplot    and    proba- 
bility   plot  are    fairly  good    for    this    set.      An    examination    of 
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the   report    indicates    that   some    bias   axists    by  the    use   of 
"average"   times    for   some   repair    time    segments   for    much   of 
the    data   instead  of   normal    data.      Tharefore,    a   further    exa   - 
mination   of  this   data    would    be   in    order. 
g.       Set    No.     15    -    AN/OYK- 14  (C) 

Set    No.    15a   -   Revised   AN/UYK-14  (C) 
Set    No.    15   illustrates   what    happens    when    the   da- 
ta  points   are  not    selected    at  random    and    thus   introduce    bi   - 
as.       In   this  case,    33    tests    were   selscted   according  to   the 
sample   selection   procedure    in   MIL-STD-471A  [Ref.    23],       How- 
ever  in   order  to   have   50   data  points    according  to    the    test 
plan    17    additional    data   points    using    "average"    repair   times 
for    several  components  from   a    previous    test    were    used.      This 
therefore    introduces    a  bias    in   the   data,    resulting   in   the 
initial   rejection    of   the    lognormal   model    (Table    3).      After 
taking   out    the    biased   data    points   the    analysis   shows    that 
the    lcgnormal    model   is  a    good  descriptor.      The  histogram 
does    not   look   very    good,    but   the    probability   does    for   the 
lognormal   model. 

This   example   then   shows    that    if    data    does    not 
fit    an   assumed    distribution    and    if   the    reason    is    bias,    then 
after    finding   a    good    reason    for    throwing   out    some    bias      po  - 
ints,    the   data    freguently   will   agree    with   a    lognormal    model. 
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h.      Sets   No   2 1    and    34   -    Fuselage    Structure 

These   two    casss    are    example   of    large    sample    si  • 
zes     (157   and    171).      For   both   data   sets    all    distributional 
assumptions   were   rejected    by   all   test    methods.      Bur   for   the 
lognormal   case    probability    plots,    boxplots    and  histograms 
are    fairly    good.      These   cases   should    be    further   invesx.iga   - 
ted. 

B.       ERRORS    IN    CALCULATED    PARAMETERS 

1 .       Error   in   MTTR    when    Assuming    an    Exponential 

Distribution    Ins tea d   of   Lognormal 

The   percent   error   in   the    mean-time-to-rspair     (MTTR) 
is  calculated   as    follows, 

IM         -    M  I 

log        exp    I 

E    =      x    100  where, 

B 

[1*0.5(7* 

M  =   the    lognormal   mean   =    2 

log 

M  =   the   exponential  mean    =    x 

exp 

The   results    of  the    percentage    srror   in   the   mean 
which    are   summarized    in   Table  5.      All    cases   havs    an   error 
less   than    18    %-      Therefore    since  the   inherent   availability 
formula   normally   used    is, 


a 

MTBF 

1 

I 

BTBF  +  MTTR 

1 

+ 

MTTR 
MTBF 
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Table    5 

PERCENTAGE    ERROR    IN     MTTB    WHEN    ASSUMING    AN 
EXPONENTIAL    INSTEAD     OF    LCGNORMAL    DISTRIBUTION 


Set 

Sample 

Sample 

Error     {%) 

No 

Mean 

Mean 

(EXP) 

(LOG) 

I3-2! 

-1- 

-2- 

-3- 

3 

1 

119.51 

112.44 

6.29 

2 

24.0 

22.  1  7 

8.25 

3 

121.63 

140.84 

13.63 

5 

57.  21 

48.6  9 

17.49 

9 

7.23 

7.22 

0.  16 

10 

9.35 

9.34 

0.08 

11 

43.05 

43.59 

1.23 

12 

15.75 

15.74 

0.02 

14 

15.73 

15.54 

1.25 

15 

15.  35 

15.59 

1.65 

17 

17.42 

17.76 

1.93 

18 

19.7 

19.1  3 

2.96 

22 

1.87 

2.0  7 

9.73 

23 

1.55 

1.4  8 

4.9 

25 

1.  76 

1.78 

1.36 

27 

0.57 

0.3  1 

1.36 

28 

0.41 

0.28 

1.29 

29 

2.45 

2.44 

0.69 

30 

1.7 

1.72 

1.45 

32 

1.64 

1.66 

0.85 

35 

4.03 

4.5  7 

11.84 

36 

2.55 

2.39 

6.96 

38 

2.26 

2.3  3 

2.96 

39 

4.  22 

4.  1 

2.73 

40 

2.5 

2.57 

2.57 

41 

3.5 

3.5  1 

0.27 

42 

3.63 

3.5  6 

0.63 

43 

2.4 

2.44 

1.78 

44 

1.72 

1.87 

7.88 

45 

2.  58 

2.60 

0.58 

NOTE:      O11I7    cases    included   when    loanormal   model   was 
accepted    (   at    least    one    out    of    three    test   methods    ) 
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and  since  for  a  reasonable  availability  MTTR  <<  HTBF  (e.g. 
MTTR  is  measured  in  hours,  HTBF  is  100fs  of  hours),  then  it 
is  shown  that  the  exponential  assumption  for  MTIR  introduces 
negligible  error  in  availability. 

C.   CONCLUSIONS 

From  the  data  analysis  in  this  study  as  well  as  the  pre- 
vious study  [Ref.  1],  it  is  concluded  that  the  lognormal 
model  is  a  good  descriptor  for  repair  times.   28  out  of  46 
seis  show  that,  within  a  0.05  level  of  significance  for  at 
least  one  of  the  3  test  methods  ,  this  model  can  not  be  re- 
jected.  Similarly,  the  data  analysis  shows  that  the  assump- 
tion of  the  exponential  distribution  should  be  rejected  in 
40  sets.   The  gamma  assumption  is  not  rejected  in  27  cases. 

The  percent  difference  in  the  MTTR,  when  assuming  expo- 
nential distribution  when  the  true  distribution  is  lognormal 
has  been  found  to  be  negligible  for  oalculated  system  avai- 
lability. 

When  the  results  of  the  statistioal  tests  indicate  oppo- 
site conclusions,  probability  plots,  boxplots  and  histograms 
are  helpful  in  determining  an  acceptable  model  and  to  explo- 
re anomalies  in  the  data.  Probability  plots  are  also  useful 
in  estimating  percentiles,  as  well  as  density  parameters. 
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APPENDIX    A 

SET    NO    1    -    FRENCH    FESSSNHEIM    PUMPS     (DOWN    TIME) 
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SET    NO    2    -    FRENCH    FESSENHEIM    PUMPS     (REPAIR    TIME) 


SAMPLE    SIZE       N    =       43 
SAMPLE    MEAN    =      24.  00 
EXPONENTIAL 
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STANDARD  DEV  =   31.95 
LOGNORMAL      DIFFERENCE 
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SET    NO    U    -    HEATER    DRAIN    PHMP 


SAMPLE    SIZE      N    =       43 
SAMPLE    MEAN    =      67.  60 
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SET    NO    5    -    MAIN     FEEDWATEF    PUMP 


SAMPLE    SIZE      N    =       33 
SAMPLE    MEAN    =      57.  21 
EXPONENTIAL 
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SET    NO    6    -    REACTOR    COOLANT    PUMP 


SAMPLE    SIZE      N    =       21  NO.     OF    CELLS      K    =       4 

SAMPLE    MEAN    =      53.57  STANDARD    DSV    =       39.22 

EXPONENTIAL  LOGNORMAL  DIFFERENCE 
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SET    NO    9    -    AN/ARC-162  (V)      (INTERMEDIATE    LEVEL) 

SAMPLE    SIZE      N    =       49  NO.     OF    CELLS      K    *    10 

SAMPLE    MEAN    =         7.23  STANDARD    DEV    =         3.44 

EXPONENTIAL  L03NORMAL  DIFFERENCE 
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SET    NO     10    -    AN/ARC-163     (7)      (INTERMEDIATE    LEVEL) 

SAMPLE    SIZE      N    =       50  NO.     OF    CELLS       K    =    10 

SAMPLE    MEAN    =         9.35  STANDARD    DEV    =         4.27 

EXPONENTIAL  L03NORMAL  DIFFERENCE 
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SET    NO    15    -    AN/UYK-14(C) 


SAMPLE    SIZE      N    =       50 
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SET    NO    15A   -    REVISED    AN/IJ YK- 1  4  (C) 


SAMPLE    SIZE      N    =       33 
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SET    NO    18    -    AN/GXS-2  (V) 


SAMPLE    SIZE       N    =       26 
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SET    NO    21     -    FUSELAGE    STRUCTURE 


SAMPLE    SIZE      N    =    1 57 
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EXPONENTIAL 


NO.     OF    CELLS       K    =       24 
STANDARD    DEV    =  1.83 

LOGNORMAL  DIFFERENCE 


PARAM1 

0.6  5 

-0.06 

PARAM2 

0.97 

MTTR 

1.54 

1.  53 

1.01 

,0 

50-TH    PERCNT 

1.07 

0.94 

13.52 

% 

90 -TH    PERCNT 

3.56 

3.33 

6.93 

% 

95-TH    PERCNT 

4.63 

4.75 

2.64 

% 

CHI-SQR    STAT 

169.06 

155.0 

DEG    OF    FREED 

22 

21 

SIGNIF    LEVEL 

0.0 

0.0 

BOXPLOT    SET21 


—  I    +         I— 


BOXPLOT    LNSET21 


**       **  OOO    0    0 


♦  I- 


00 
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HISTOGRAM    SET    21 

EACH    *    REPRESENTS 


2    OBSERVATIONS 


MIDDLE  OF 

NUMBER 

OF 

INTERVAL 

OBSERVATIONS 

0. 

33 

***************** 

1. 

73 

*****  **  **  ******  ************  *****  ***** 

2. 

27 

*****  ********* 

3. 

7 

**** 

4. 

4 

** 

5. 

2 

* 

6. 

5 

*** 

7. 

3 

** 

8. 

1 

* 

9. 

0 

10. 

2 

* 

TOGRAM  LNSET21 

MIDDLE  OF 

NUMBER 

OF 

INTERVAL 

OBSERVATIONS 

-2.5 

3 

*** 

-2.0 

2 

** 

-1.5 

12 

*****  ******* 

-1.0 

16 

**************** 

-0.5 

34 

********************************** 

0.0 

36 

************************************ 

0.5 

27 

*************************** 

1.0 

10 

*****  ***** 

1.5 

7 

*****  ** 

2.0 

8 

*****  *** 

2.5 

2 

** 
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SET    NO    28    -    MAIN    ROTOR  HEAD    AND    BLADES     (ELAPSE    TIME) 

SAMPLE    SIZE      N    =       33  NO.     OF    CELLS       K    =       6 

SAMPLE    MEAN    =         2.07  STANDARD    DEV    =  1.87 

EXPONENTIAL  LOGNORMAL  DIFFERENCE 


PARAM1 

0.4  8 

0.46 

PARAM2 

0.52 

MTTR 

2.07 

2.05 

0.89 

50-TH    PERCNT 

1.1*4 

1.58 

9.09 

5? 

90 -TH    PERCNT 

4.77 

4.00 

19.31 

n 

95-TH    PERCNT 

6.21 

5.20 

19.33 

r 

CHI-SQR    STAT 

11.91 

2.09 

DEG    OF    FREED 

4 

3 

SIGNIF    LEVEL 

0.  180  E- 

■01 

0.554E+ 

00 

BOXPLOT    SET28 

—  I      «■  I- 

BOXPLOT    LNSET28 
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HISTOGRAM    SET28 


MIDDLE    OF 

NUMBER 

OF 

INTERVAL 

OBSERVATIONS 

0. 

0 

1. 

16 

** ************** 

2. 

6 

*****  * 

3. 

7 

^ft  Jf<    -^S.  ^*   t*     ^l   »^* 

4. 

2 

** 

5. 

1 

* 

6. 

0 

7. 

0 

8. 

0 

9. 

0 

10. 

0 

11. 

1 

* 

TOG? AM    LNSET2  8 

MIDDLE    OF 

NUMBER 

OF 

INTERVAL 

OBSERVATIONS 

-0.8 

3 

*** 

-0.4 

2 

** 

0.0 

8 

*****  *** 

0.4 

6 

*****  * 

0.8 

6 

*****  * 

1.2 

6 

*****  * 

1.6 

1 

* 

2.0 

0 

2.4 

1 

* 
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SET  NO  34  -  FUSELAGE  STRUCTURE 


SAMPLE  SIZE   N  =  171      NO.  OF  CELLS   K  ■  34 
SAMPLE  MEAN  =    2.83      STANDARD  DEV  =    4.60 

EXPONENTIAL     LOSNORMAL      DIFFERENCE 


PARAM1 

0.35 

0.36 

PARAM2 

1.26 

HTTR 

2.8  3 

2.68 

5.26 

% 

50-TH    PERCNT 

1.96 

1.43 

36.77 

or 

90-TH    PERCNT 

6.5  1 

6.03 

7.95 

.'O 

95-TH    PERCNT 

8.46 

9.05 

6.51 

% 

CHI-SQR    STAT 

126.  65 

92.85 

DEG    OF    FREED 

32 

31 

SIGNIF    LEVEL 

0.596E- 

•07 

0.596E-07 

OT    SET34 

1+       I *    **O0  O    300 

30XPL0T    LNSET34 


** 
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HISTOGRAM    SET34 

5    OBSERVATIONS    ARE    BELOW    THE    FIRST    CLASS 

MIDDLE    OF  NUMBER    OF 

INTERVAL  OBSERVATIONS 

0.500  38  ************************************** 

ll 000  33  ********************************* 

lis  00  25  ************************* 

2.000  15  *************** 

2.500  9  ********* 

3.000  13  ************* 

3-500  2  ** 

U.  000  3  *** 

4.500  3  *** 

5.000  3  *** 

5.500  0 

6.000  5  ***** 

6.500  1  * 

7.000  1  * 

HISTOGRAM    LNSET34 

MIDDLE    OF  NDMBER    OF 

INTERVAL  OBSERVATIONS 

-2.500  2  ** 

-2.250  0 

-2.000  0 

-1.750  3  *** 

-1.500  5  ***** 

-1.2  50  0 

-1.000  11  *********** 

-0.750  15  *************** 

-0.500  6  ****** 

-0.250  13  ************* 

0.0  19  ******************* 

0*250  15  *************** 

0.500  19  ******************* 

0^750  12  ************ 

1.000  18  ****************** 

1.250  3  *** 

1. 500  8  ***** *** 

1.750  6  ****** 

2.000  2  ** 

2.  ">5Q  7  ***** ** 

2.500  1  * 

2.750  3  *** 

3.000  0 

3.250  2  ** 

3.500  1  * 
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SET    NO    41     -    MAIN     ROTOR  HEAD    AND    BLADES     (MAN    HOURS) 

SAMPLE    SIZE      N    =       35  NO.     OF    CELLS      K    =       7 

SAMPLE    MEAN    =         3.50  STANDARD    DSV    =         3.42 

EXPONENTIAL  LOSNORMAL  DIFFERENCE 


PARAM1 

.0.29 

0.90 

PARAM2 

0.71 

MTTR 

3.50 

3.51 

0.27 

% 

50-TH    PERCNT 

2.43 

2.47 

1.57 

% 

90-TH    PERCNT 

8.06 

7.25 

1  1.29 

% 

95-TH    PERCNT 

10.49 

9.83 

6.72 

% 

CHI-SQR    STAT 

11.20 

7.20 

DEG    OF    FREED 

5 

4 

SIGNIF    LEVEL 

0.  476 E- 

01 

0. 126E+ 

00 

BOXPLOT    SET41 

-I       +  I 

3OXPL0T    LNSET41 
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HISTOGRAM    SET41 


MIDDLE    OF 

NUMBER 

OF 

INTERVAL 

OBSERVATIONS 

0. 

5 

***** 

2. 

13 

*****  ******** 

4. 

10 

********** 

6. 

4 

**** 

8. 

1 

* 

10. 

1 

* 

12. 

0 

14. 

0 

16. 

0 

18. 

1 

* 

TOGRAM    LNSET41 

MIDDLE    OF 

NUMBER 

OF 

INTERVAL 

OBSERVATIONS 

-0.5 

3 

*** 

0.0 

7 

*****  ** 

0.5 

7 

*****  ** 

1.0 

3 

*** 

1.5 

9 

*****  **** 

2.0 

5 

***** 

2.5 

0 

3.0 

1 

* 
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APPENDIX  B 
Table  6 
TOTAL  RESULTS  (46  cases) 
Exponential  Distribution  :   Only  Chi  Square  test  used 

Accept      Rejert 
|    6      |    40     | 


Lcgnormal  Distribution  :   Chi  Square,  Kolmogorov  Smirnov 

and  W  test  used. 

K^S  -est 


-• 

W    i 

A 

A 

It 

R 

A 

R 

te 

st 
R 

No    H 

K-S 

A 

tes- 

test 
R 

■ 

i 

A 

R 

1 
1 

12 

1 

7 

1  I 
i 

Ch. 

Sg. 

!St 

A| 
I 

a 

Chi    Sq. 

2 

0 

I       2    | 

test 

1 
I 
1 

1 

2 

3 

7| 

I 

I       8    | 

(  34  cas  es  )  (12  cases  ) 

Gamma  Distribution  :   Kolmogorov  Smirnov  test  used. 

M.O.H.  estimators 
A  R 


M.L.E.      A   |  19    |  3 

L       .1.  _i 

i         i  i 

estimators   R       J           0      j  19      | 
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Table    7 


PUMP    DATA    RESULTS     (7    S2tS) 


Exponential    Distribution    :      Only   Chi    Square   test    used 

Accept  Reject 


I  1 


Lcgnormal   Distribution    :      Chi   Square,    Kolmogorov    Smirnov 

and   W   cast   used. 


K-S    test 


A 
R 

-• 

W    1 

A 

A 

4- 

R 

W 
A 

R 

test 

R 

tes*£ 

A 
R 

No    W 

K-S 

A 

test 

test 
R 

! 
I 

I 

! 
I 
I 

0 

1 

2 

I 

0    I 
i 

Chi    Sq. 

0 
0 

1       °    1 

test 

0 

1 

0 

2! 

I 

1 1 

1       1    I 

(    6    cases  )  (    1    case    ) 

Gamma   Distribution   :      Kolmogorov   Smirnov   test    used. 


B.L.E. 


M.O-M.  estimators 
A  R 

A   I     0       1 


estimators  R 


0   |    6 
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Table    8 
REPORT    DATA    RESULTS     (13    sets) 
Exponential    Distribution    :      Only    Chi    Square   test    used 

Accept  Reject 


Lognoraal   Distribution    :      Chi   Square,    Kolmogorov    Smirnov 

and    H    '-.est    used. 


k-s   test 


Chi    Sq 
test 


a 


A 

R 

H    t 

es 

4- 

W 

te 

-.4- 

A 

R 

A 

R 

i 

I 

4 

0 

2 

o  ! 

I 

! 

! 

0 

0 

2 

3 

Ch. Sq. 

test 


No    W    test 


K-S    test 
A  R 


Al        0        |        1     I 

I 1 i 

R!       0       |       1    1 


{    11    cases   )  (    2    cases   ) 

Gamma   Distribution   :      Kolmogorov   Smirnov   test    used. 

a.o.a.    estimators 

A  R 


B.L.B. 


estimators   R  3  3 
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Z§k2£   1 
FIELD    DATA    RESULTS     (26    sets) 

Exponential    Distribution    :      Only   Chi    Square   test    used 

Accept  Reject 

I 

22  J 


Lognormal   Distribution    :      Chi   Square,    Kolmogcrov    Smirnov 

and    »    -est    used. 


No    W    test 


Chi    Sq, 
test 


8 


K- 

•s 

test 

H 

A 

A 
te: 

3t 

R 

R 
W    test 

A         a 

! 

\    3 
1 

0 

3            ,    I 

I 

1    1 

1 

1            2| 

Ch. Sq. 
tesr 


K-s   test 

A  R 


A  I  1 

R|       0  6 


(    17    cases    )  (    9    cases   ) 

Gamma   Distribution   :      Kolmogorov   Smirncv   test    used. 


H.L.E. 


M.O.M.  estimator: 
A  R 

A  11  |  5 

I  I 


estimators   R 


0       1        10 
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Table  10 


TOTAL  %    RESULTS  (w.r.t.  sample  sizes) 


For  exponential  distribution  : 
Sample  Size         A 


20-29 

30.8 

59.2 

30-39 

8.3 

91.7 

40-50 

0. 

100. 

GT.  50 

8.3 

91.7 

For  lognormal  distribution 


Sample 
Size 


20-29 
30-39 
40-50 
GT.50 


Chi- 
A 

Sq 

.test 
R 

W 
A 

test 
R 

K-S 

A 

B  " 

53.3 

46.2 

71.4 

28.6 

46.2 

53.8 

66.7 

33.3 

58.3 

41.7 

50. 

50. 

58.3 

41.7 

58.3 

41.7 

33.3 

66.7 

25- 

75. 





16.7 

83.3 

For   gamma   distribution    : 


Sample 

Size 

A  in 

both 

84.6 

R  in 
both 

7.7 

A  in  HLE 
R  in  MOM 

20-29 

7.7 

30-39 

33.3 

5  0. 

16.7 

40-50 

25. 

52.5 

12.5 

GT.50 

16.7 

15.7 

66.6 
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APPENDIX  C 

COMPUTER  PROGRAMS  FOR  DATA  ANALYSIS 

1.  LGNRML  (FORTRAN)  [Ref.  1] 

The    program    makes    a   Chi-Sguared   Goodness-Of-Fit    test    for 
the    exponential   and   lognormal   assumptions.      The   main    program 
computes  expected   values    of    ordered   observations,    sample 
mean    and   standard   deviation    based   on    tha    maximum    likelihood 
estimates   of   the   parameters    of   the    lognormal    distribution, 
calculates    the    50thr    90th,     95th    percentiles   of   the   exponen- 
tial   and   lognormal    distribution    assumptions,    percentage 
differences   between   the   calculated    parameters   of    the 
exponential  and    lognormal   distributions. 

2.  BOXPLOT     (MINITA3)     [Ref.     5] 
Syntax:      Boxplot    Vector 

Description:  Generates  a  boxplot  display  for  a  vector 
of  data.  A  rectangular  box  with  ends  corresponding  to  lower 
and  upper  quartiles  is  presented  with  the  median  marked  with 
an  plus  sign,  and  lines  are  drawn  on  each  side  of  the  box 
with  crosses  marking  the  lowest  and  highest  data  values  wit- 
hin half  guartile  distance  in  each  side.  Outliers  are  mark- 
ed   with    stars   and   circles. 

3.  NORMPLOT     (APL)     [Ref.     21  ] 
Syntax:       Normplot    vector 

Parameters:    Wid-   controls   the   horizontal    size    of    display 
Dep-    controls   the    vertical    size    of   display 

Subprograms:      NSCORSS,    PLOT 

Description:      Generates    a   normal    probability    plot    for    a 
vector   of    data.    If    data    fits   normal   then    plot    is    straight 
line. 
U.       SXPONPLOT     (APL)     [Ref.    21] 

Syntax:      Sxponplot   vector 

Parameters:    Wid-    controls  the   horizontal   size    of    iispiay 
Dep-    controls   the    vertical    size    of    display 

Subprogram:       PLOT 
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Description:      Generates    a   exponential   probability   plot 
for    a    vector   of    data.      If   data    fits    exponential   then    plot 
is   straight   line. 

5.  GAM     (APL) 

Syntax:      GAM   vector 

Subprograms:       ESTGAM,    GAMMA 

Description:  Perform  K olmogorov-Smirnov  Goodness-of-Fit 
test  and  gives  test  statistic  and  "reject"  /  "do  not  reject" 
idea    for   gamma   assumption   for  both   MOM   and    MLE   estimators. 

6.  KOL     (APL) 

Syntax:      Kol   vector 

Subprogram:       MOM,    INVV 

Description:      Performs    a   Kolmogorov    Smirnov   Goodnees-Of- 
Fit    test    for   a    data    vector    for    lognormal    assumption   and 
print   test    statistic.      If   test   statistic   is   bigger   than  the 
tabulated   value   "reject."    the   assumed    distribution,    otherwise 
"do    not   reject". 

7.  W    (APL) 

Syntax:   A  9   Vector 

Description:   Performs  a  W  test  for  the  lognormal  as  - 
sumption.  A   must  be  a  data  vector  (from  Table  IX,  Ref.  2) 
and  X  vector  that  is  testing  against  lognormality .   A  and  X 
must  have  the  same  dimension,  output  will  be  W  test  statis- 
tic.  Then  using  the  value  of  A  ,/u  and  €  (Table  XI,  Ref.  2) 
for  the  appropriate  sample  size  and  than  using  standardized 
normal  distribution  determine  the  probability  of  obtaining 
value  less  than  or  eg ual  Z,    which  is  the  significance  level 
of  the  test. 
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cm 

C23 

C3: 

L'43 

cs: 

C7J 

c?: 

ClOJ 

till 

L121 
C133 
C143 
CloU 
1161 
C17J 
E13D 
C193 
C203 
C213 
C222 
C23D 
C243 
C253 
C261 
C273 
C28D 
C293 
C30  3 
C313 
C323 
C333 
C343 
C353 
E36H 
[373 
C383 


KOLMOGOROV  5MIRNOF  TEST  FOR 
FOP:  MOM  AW  It  tfLE  ESTIMATORS 


v  gam  can? 

7  GAM  V 
A  PROGRAM  PERFORMS 
ft   GAMMA  ASSUMPTIOW 

ESTGAM  Y 

xi«.Y^-eciJ2D 

>(2<-i^ec2  J  2: 

N  <-  f  Y 

24-G*886*«*0*5 

«2*-EC2J  1] 

K«-Of  «*«■<!  (N-l)  )tN 
i«-ii«-0 

j<-jj<-0/0 
L2J  K-I  +  l 
■»L.lX\(I>H) 
Al   GAMMA  X1QI] 
JfJfP 

-n-2 

R  RfR[  f  R  ] 
1   TEST  STATISTIC  FOR 

•   REJECT  FOR  ALFHfi  =  0.05  FOF: 

*i-5 

L3J'   DO  WOT  REJECT  FOR  AL.PHAxO»05  F  °  p: 
l_5t  IIfII  +  i 

-tL^X  \ < II >N) 

A2   GAMMA   X2C11! 

JJf JJ,P 


MOM  ESTIMATORS   IS   '  »  *•  R  ft  C  1  3 
MOM  ESTIMATORS   ' 

40M  ESTIMATORS 


>*L! 


K) 


L^JS«.|  (  JJ 

S5f,5[  f  S3 
1   TEST   STATISTIC   FOR   ML.E  ESTIMATOR 

41-7x1  (ssci3  (Z) 

'       REJECT      FOR      AL.FHAxO.05      F  ° R      ML-E      ESTIMATORS 
■499 
L-7J   ' 

->99 


t55[l 1 


DO  WOT  REJECT  FOR  A  l_  P  H  A  x  Q  ♦  0  5  F  °  p:  M  L  E  ESTIMATOR 


GAM  5ET1 
TEST  STATISTIC  FOR  MOM  ESTIMATORS   IS  0,363 
REJECT   FOR   Al_PHA=Of05   FOR   MOM   ESTIMATORS 
TEST  STATISTIC  FOR  MLE  ESTIMATOR   IS  0*154 
REJECT  FOR  Al_PHAxO»05  FOP'  M  L  £  ESTIMATORS 
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[13 

C23 
[33 

C43 

C53 

161 
Z71 
LSI 

Zfl 

C103 
C113 
L123 
C133 
C143 
C153 
C163 
[173 
C183 
C193 
C203 
C213 
[223 
[233 
C243 
C253 
C263 
[273 
[283 
[293 
C303 
C313 
[323 
[333 


v    kol  CD lv 

7   KOL   X 
A  PROGRAM  PERFORMS   K OLMOGOR O V - 5 M I R M O V 
fl  GOODHE55-OF-FIT  TEST  FOR  LOGI-IORMflL 
ft  ASSUMPTION   ,   X  MUST  BE  DATA  ARRAY 
#      ORDER  DATA  AND  FIND  SAMPLE  MEAN  AND 

x«.xL*x3 

MOMfX 
ft   PUT   DATA   IN   STANDARD   NORMAL  CASE 

2<-(  (  »X)-XB  )  iS 
NffK 

ppp«-0/Q 

ifii<-0 
L3; iifi i+i 

4  L  4  X  l  (  I  I  >  N  ) 
ft  FIND  P  VALUE  FOR  EVERY  DATA  POINTS  AND 

l-i:  K-I  +  l 

->L4X  \  (  I>??9) 
IHVV   P<-IX0,001 
-»L2X\  (2Z>2CH3) 

L2 ;ppf>pfp f  P 
■H-3 

L4:Kf(t(N-l) )rN 

K«-0>  X 

CCf (fK)-(fPPP) 

DC<-<  ,  CD) 
PPP*-PPP,DC 
H     FIND  MAXIMUM  ABSOLUTE  DEVIATION  AND  PRINT 
KK«-|<PPP-K) 

L<-KK[fKK] 
■TEST   STATISTIC  ='»tLC13 

tv«.o.886*<n*0*5) 


VARIANCE   OF   LOG   DATA 


PUT   I N  A  R  R  A  Y 


AS   TEST   STATISTIC 


TABULATED   VALUE 


•  TV 


KOL   S  E  T  1 
TEST   STATISTIC   =0,146 
TABULATED  VALUE  =0.135 
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1X1 

[23 

[33 

[43 
[53 
[63 
[73 
[83 
[93 

[ion 

Ell] 
C  123 
C133 
[143 
[153 
[163 
[173 
[183 
[193 
[203 
[213 


v  w  CD 3 ^ 

7  V  W  K J Ai  J  A2  J  A3  J  S2  J  K  J  N  J  NN 

ft   PRODUCES   W   TEST   STATISTIC   FOP:  W   TEST 
n   fsrVECTOft  FOR  W  TEST  LENGTH  OF  ARGUMENT 
A  X  =  VECTOR  THAT   IS  TESTING 
ft  X  AND  Y  MUST  BE  SAME  DIMENSIONS 

N«-f  X 

X4-&4-0 
ft  ORDER  DATA 

x«.kc*X3 

ft  CALCULATE  S*2  TERM 
A  J  4-  +  /  <t  X 

ft2<-  +  /  <•«)  *2 
A3<-<  «1  *2>  *n 
S2<-A2-«3 

ft  CALCULATE  B  TERM 
K  f.  L  H  -r  2 

L  x  <-  »  x 

LL;«.u;;[tLX] 

B  *-  +  /  T'  X  (  (  K  f  L  L  X  )  -  (  K  ^  L  X  )  ) 
ft  FIND  W  TEST  STATISTIC 

2<-(»*2)^-52 

1   W  TEST  STATISTIC  =   ' 


AND   PRINT   AS   W 


r  t' 


A43   W   SETl 
TEST  STATISTIC 


=  0,967 
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C13 

[2: 

[33 

C43 

C53 

C63 

C73 

C33 

C93 

C  1 0  3 

C113 

C123 

C133 

C143 

C153 

C163 

H173 

C183 

C193 

C203 

C213 

C223 

C233 

C243 


v 

7   A 
ft   IN 

ft  GI 
A  A 

a<-L 
<=•(-! 
-♦Ox 
R<-  ( 
v<-< 

->L  X 
T>V 

-♦(1 

-►L 
FrjV 

E  +  2 
z<-a 

-»<  1 

+o 

7 


GAMMA   [03? 

gamma  ;■; ;  V  J  Y  J  z  ;  R  J  A  ;  N 
C,   GAMMA  FMC,-5HAPE=A>0,   VECTOR  ARG,  X > Q ,   SCALE  PAftAM 
VE5  THE  LEFT  TAIL  AREA  FOR  GIVEN  A  AND~X 
IS  SHARE  PARAMETER,   X  IS  DATA  VALUE 
A 

-  *  -  X 
I  1=A 

X  <  A  )  v  ■<  i  1 

/X 

~7  +  «<  r/Y)-rA+l 

+  +  /  X  \  Y  o  ,-^A+\H 

f3  +  QLC)  xil  <A 

f\t X]fZX(Y*A)i(  ;  A)  XX Y 

/\f  K3<-<  Y-A)  X2X  (  X/  (  Yx  *-  t  -rfi-1  )  o  ,  t- A  -  \  ft- J  )  x  (  Y*A-ft  ) -r  i  A- A 

X \0=+/"v 
*V)/X 

0T2xa 

CF  Y 

+  3  +  01-C)  x  t  1  <  ft 

*V)/\fK]<-l-ZxCf»A)x(X-T)r!A-l 

*V)/lfX1fl-TxZx(X/(Yx  *-  Y^.ft-1  )  o  ,  ift-^ft-J  )  x  (Y*A-ft)  +  i  A-fi 


=  1 
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C23 
C33 
C43 
C53 
C63 
C73 
C83 


ft  ' 

ft  M 

B<- 

X«- 

«1 

X"> 


IN 
INV 
INV 
ORM 
I  (1 
(-• 

♦-2. 
«-l  + 


vv 

V  A 

V  A 
AL 

X  A  > 
BXB 

515 
xx  1 

l*x 

-2x 


CD3* 


'  R 
DI  5 

0*5 
)  tO 

517 
.43 
2>- 

ZZ  x 


ET 
TR 
)- 

♦  5 
-t-x 

27 
x 

A> 


URNS  THE  VALUE  OF  STftUDftftD 
IBUTION  AS   LEFT   TAIL   AREA   A 
A 

x0. 302353+0. 010328*- 
88+xxO. 189269+0. 001308xx 


0. 


v  mom  CD3^ 

<?  MOM   X 
C  1  3    fl  RETURNS  SAMPLE 

C23    x»«-<+/x)rfX 

C33    s<-(+/(x-xb)  *2)+' 


MEAN     *HZ>     VARIAIICE 


I+f 
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cm 

C23 
C33 
C43 
C53 
C63 
C73 
C83 

c?3 

C103 

cm 

C123 
C133 


CD3v 


<7  E3TGAM 

V  E3TGAM  X'JAJBJFJFP 
n  RETURNS  METH,   OF  MOMENTS  EST,   THEN  MAX  LIK   EST,   OF  GAMMA  DATA 
ft  OUTPUT  AS  2*2  MATP:I*   IT5  FIRST  ROW  MOM  AND  SECOND  ROW  MLE 
ft  ESTIMATORS  5HAPE  AND  SCALE  PARAMETERS  RESPECTIVELY 

MOM  A 

ftfXBrBfSTXB 

E<-A,B 

»<-<  t/tX)rfX 
L  ♦  A  Q  «-  A 

FfBf  (  iArXB  )  -0   f^1   A 

FFf (rA)-l   PSI   A 

fl(.fl-F.FP 

-»"-x\  (0*0001)  <  l«-«0 


E<- 


2  2   P  E  ,  A  ,  X  B  -  A 


v  rsi   can? 


CI] 

C23 

C33 

C43 

C53 

C63 

C73 

C83 

C93 

C103 

C113 

C123 

C133 

C143 

C153 

C163 

CI  7 3 

C  133 

C193 

C203 


P«-N 


PSI 
THE 


IViJIVJKJKK 


I  T 


ZJ  T 


ORDER  OF  THE  DERIVATIVE 


IS   THE  ARGUMENT, SCALAR   OR 


ft     N     ; 

fl      *'      (  >  0  ) 

C<-10 

I V  *.  \  f  Y  «-  ,  Y 

P«.Z«-K*.(fY)f0 

KKffC-'l'CJIVf  (VfY(C)/IV] 
4LlX\<fIV)af(#V)/XV 

T<-Y[JIV] 

I+-0 

L2 :  i<-i  +  i 

Ti>KK[I]fT[I] 

zCi3«-(  JN)x+/<  <  ff-i)  +  lKKCi3)*-l+w 

->L2X  \  I  <  P  jiv 
Z«-V\Z£lf  JIV] 
KfV\KK 

L-1J-*S1X\N>0 

P*-(«KrY)-{2xK+r ) *-l 

+2+ai-c 

51  {  P«.  (  (  l  N-l  )  x  (  T  +K  )  ft-N  )  + 
F'M  (~1  )  *»•<•+•!  )  XR  +  Z  +  N     JEX 


OF      THE 
V  E  C  T  O  R 


PSI      FUHCTIOM 


I  H  )  X  0 


X <  Y+K ) *-N+l 
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